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Abstract. We examine a family of microscopic models of plasmas, with a 
parameter a comparing the typical distance between collisions to the strength 
of the grazing collisions. These microscopic models converge in distribution, in 
the weak coupling limit, to a velocity diffusion described by the linear Landau 
equation (also known as the Fokker-Planck equation). The present work ex- 
tends and unifies previous results that handled the extremes of the parameter 
a, for the whole range (0, 1/2], by showing that clusters of overlapping obsta- 
cles are negligible in the limit. Additionally, wc study the diffusion coefficient 
of the Landau equation and show it to be independent of the parameter. 

1. Introduction 

Particles in a plasma experience grazing collisions because they are ionized, in- 
teracting even at long distances as described by the Coulomb potential. So far, the 
full Coulomb model has been impossible to handle rigorously in a scaling limit (see 
[L] for a heuristic argument, and [DP] for a rigorous partial result), and the strategy 
has been to use an approximation by soft-sphere models with their bump-function 
potentials. 

Microscopically, soft-sphere models consist of a lightweight particle traveling 
through a random configuration of large stationary particles, called "obstacles" 
or "scatterers," whose shape and density are determined by a parameter. The 
lightweight particle grazes obstacles when it gets within their ranges of influence, 
called "protection" disks. In two dimensions, this can be visualized as a ball rolling 
through a random field of hills. It is desirable to understand these microscopic 
models in the weak coupling limit, when the radius of the obstacles goes to zero, 
and to derive rigorously a macroscopic description in terms of a linear PDE. 

More precisely, we introduce a parameter a € (0,1/2]. The i-th obstacle, centered 
at ri (see figure 1), is described by a suitably smooth and compactly supported (say, 
in the unit ball) radial potential V whose rescaling is: 

V^ix~n) = e"V{\x^n\/s). 

We use the convention of [DR]; another convention for rescaling has e^, as in [DGL]. 
Then the microscopic dynamics for an obstacle configuration, w, are Newtonian: 

[ = <^ <(0) = X, 

(1-1) \ «s = -E^^o«-^')' <{^) = v. 
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Figure 1. Example of the soft-sphere model: a configuration of 
obstacles and the corresponding trajectory of the light particle. 



The collection of obstacle centers, uj := {ri : i E Z}, is a realization of the Poisson 
point process in with intensity ;= e~'^°'~^p, e.g., the expected number of 
obstacles in A c is E'^(iV(A)) = p\A\s-'^"-'^ . 

The case a = 0, where the particle travels a relatively long distance between 
collisions and each obstacle has a relatively large influence, corresponds to the 
Boltzmann-Grad limit of the hard-sphere model of a dilute (or Lorenz) gas, with 
the macroscopic evolution given by the Boltzmann equation. (See, for instance, [SI] 
and [BBS]; and for the quantum Lorenz gas, [EYl].) 

For a G (0, 1/2], the number of obstacles must be about £-20-^+1 p^j- ^nit of 
volume in order to have a net effect that is nonzero and finite (this is more than 
in the Boltzmann-Grad limit, where the number is about £^~'^). This large number 
of obstacles balances the small factor of in the rescaled potential, and as e ^ 
in this, the weak coupling limit, the macroscopic dynamics of these plasma models 
are given by the linear Landau equation. 

Previously, the results for this family of models were incomplete. Desvillettes 
and Ricci proved a weak version (convergence in expectation) of the two-dimensional 
weak coupling limit for < a < 1/8 [DR]. They approximated the Landau equation 
by the Boltzmann equation in order to use a modification of Gallavotti's technique 
for the Boltzmann-Grad limit. From the outset, however, they needed a to be 
small, so that the radius of each obstacle is much smaller than the expected free 
flight time. The method is further limited to the regime < a < 1/8 by an estimate 
of the probability of self-intersection, but can be improved to include < a < 1/4 
(see Appendix). 

On the other hand, Kesten and Papanicolaou proved a stronger convergence 
(in law) of the weak coupling limit at the upper endpoint of the parameter range, 
a = 1/2, for dimensions three and higher and with a general random field, called the 
stochastic acceleration problem [1\J'1]. Then Diirr, Goldstein, and Lebowitz proved 
a two-dimensional version, with the Poisson distribution of obstacles [DGL]. More 
recently, Komorowski and Ryzhik handled the stochastic acceleration problem in 
two dimensions [KR2]. (See also [EY2] and [IlLW] for the quantum weak coupling 
limit.) 
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The difficulty in the middle range of a is that obstacles overlap a great deal more 
than for small a, but their individual influence is greater than for a = 1/2. There's 
an additional difficulty in the two-dimensional case because the probability of self- 
intersections is nontrivial. Fortunately, it turns out that "bad" self-intersections 
(ones that are repeated or almost tangential) are negligible in the scaling limit. 
This is important because such self-intersections cause a correlation between the 
past and the present (see Figures 2 and 3), and controlling memory effects is the 
main difficulty in these problems. In higher dimensions, this is unnecessary, because 
the probability of any self- intersection is negligible in the limit. 




regime of 
dependence 



Figure 2. Correlation within a trajectory due to a small-angle 
self- intersection. 



many self-intersections 




regime of 
dependence 



Figure 3. Correlation within a trajectory due to many self-intersections. 



Theorem 1.1. For a G (0, 1/2), the family of stochastic processes {v^{t))t>o con- 
verges as £ ^ Q in law to the velocity diffusion (v{t))t>o generated by A-u, the 
Laplace- Beltrami operator on := {w E : \w\ = |wo|}- 

In particular, z//o(x,w) is an initial distribution of positions and velocities, and 
ow for the microscopic dynamics (1.1), then 





'^a,u,e «S the 



(1.2) 



f^it, X, v) = nfoi^'i^, v))] ^ hit, X, v), 



where h is the solution of the linear Landau equation: 

j {dt-hv- S/^)h{t, X, v) = CAi,/i(t, X, v); 
I h{0,x,v) = fo{x,v). 



(1.3) 
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Moreover, the microscopic distinctions of the obstacles' steepness and density all 
disappear in the scaling limit, and the models have the same macroscopic behavior: 

Proposition 1. The diffusion coefficient in (1.3), C, is independent of a G (0,1/2] 
and can be expressed by the following formula: 

The main novelties in the present paper are better estimates on clusters of over- 
lapping obstacles and the amount of time spent interacting with each cluster. By 
controlling these quantities, we show that the total influence of clusters is negligible 
in the scaling limit. The outline of the paper is as follows: 

• To prove convergence in law, we must show that the measures induced 
by the microscopic processes {v^{t)) on C([0, T]^; x R^) converge weakly 
to I' = z^Q induced by the diffusion {v{t)). 

— To this end, we first define some stopping times to eliminate wild be- 
havior (e.g., too many self-crossings, or at too small an angle), which 
is negligible in the limit (Section 2). 

— Next wc prove that the family of stopped processes is tight (Section 
3). 

— And we identify the limit as a velocity diffusion (Section 4). 

• Then we examine, in Section 5, the particular case of convergence in ex- 
pectation, using a Gallavotti-type method and a combinatorial argument 
about clusters of obstacles. 

• Finally, in Section 6, we prove the Proposition about the diffusion coefficient 
being constant in a, and show that the formulas of [ : .'] and [DGL] agree. 

• The Appendix contains a narrower extension of the techniques of [DR]. 

2. The cut-offs 

We prove the theorem for the processes with stopping times that can be removed 
for the limiting process. These cut-offs prevent wild behavior that would lead to 
a correlation (the regimes of dependence in figures 2 and 3) between the past of a 
trajectory and its future. First, for p € C([0, r]^; x R'^)^ define: 

(2.1) Q{t) := f p{u)du. 

Jq 

Then we define the stopping times as follows. The first time the trajectory ap- 
proaches its past within distance a and with angle less than is cut off by r^^a- 

T^^a := inf{t :3s G [0,<] such that \Q{s) - Q{t)\ < a, 

^^■^^ min »(s) -piu) < 0, and /^""jl > cos0}. 

«e[s.t] \p{t)\\p{s)\ 

The first time trajectory crosses itself more than K times is cut off by tk- 

Tx := infjt > :3sj < i,;, for i ~ \, . . . and ti < ■ ■ ■ < tj^ = t, 

^ ■ ^ such that Q{si) = Q{ti)\fi} 

Additionally there is a velocity cut-off, to prevent the trajectory getting "stuck" 
somewhere or going too fast: 

(2.4) T.„ := inf{t > : \\pit)\ - vo\ > vo/2}. 
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We call the overall stopping time r: 

(2.5) r := niin{r0, a, Tif,T„}. 

These stopping times actually depend on e, so we write the stopped processes as: 

(2.6) :=<(iAr^). 

And their induced measures are written as z^^. To see that the cut-offs can be 
removed, namely that v"^ and the measure induced by the original process, v"^, have 
the same limit, v, we claim that, for arbitrary T > 0: 

(2.7) limi.^(T<r) = 0; 

(2.8) lim z/(T0,a < T) = 0; 

0,a— i-0 

(2.9) lim v{tk <T) = 0. 

K—Kx 

The first statement is immediate. The second follows from a self-crossing lemma 
in [DGL, pp. 228-9], which says roughly that tangential self-intersections are neg- 
ligible. The third follows as a corollary, because the trajectory is almost surely a 
continuously differcntiable curve. 

3. Tightness 

We will show that the family of processes is tight, to get the weak convergence 
of the induced measures i>^. First we need some generalizations of the lemmas in 
[DGL, § 4]. 

3.1. The Martingale-compensator decomposition. We can view the Poisson 
point process of obstacles from the particle's perspective (a locally Poisson process 
denoted N^, realizations of which are measures on 5* x M+), and the difference 
(denoted M^) turns out to be a martingale. This is expressed in the following 
equation, where the first term on the right-hand side is a martingale, and the 
second is a measurable left-continuous process called the compensator. 

(3.1) 

' N'^{da,du)f{a,u)^ f f M'Ma,du)f{<j,u)+[ f p%{da,du)f{a,u). 

JSJO J SJQ 

Then we write the potential as F^{y) := — X]^^a(2/)' the stopped position 
process x^'" = x^(t A r^), and a;^'_^(t) := — x^^ + ae. It is convenient to express 
Wj'" in terms of N^,: 

(3.2) vr-vo= / K{da,du) F^«;:(s))rfs + Ai;r(0- 

JR2 Jo Ju 

3.2. Bounding clusters of overlapping obstacles. To control the probability 
of A scatterers being found within distance ne of the particle at time t, we introduce 
some notation: 

Slix) :=i?(x,ne), 
AAr(a:) := |c^n5^(x)|, 
B := sup|y(a;)|. 

Lemma 3.1. For any A,e,T > 0, 

P^({supAA,^'"(xr) > A}) < f ^ + ^^^"-^^ + g32„^pe-— -A_ 

t<T \ ne J 
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Proof. The proof follows [DGL] with appropriate modifications. We start from a 
result of the conservation of energy, with vq = 

< 1 + 4e"supAAf'"(a;^'")B. 

s<t 

Hence < 1 + 4e"AB, if A is a bound on AAf'"(x|'"). 
Then, if we set tx := inf{t > : Af,'^'" {x^'") > A}, we have 

{ta > T} c {sup<'" < 1 + 4\e"B} c {supx^'" < (1 + 4Ae"B)r}. 

t<T t<T 

That is to say, if the particle does not come close to A scatterers at once before time 
T, then < l + 4Ae"B for all t before T; hence x^'" is bounded by (l + 4Ae"B)r. 

Now we consider the square [-(1 + 4Ae"B)r, (1 + AXe°'B)T]'^, which we call T. 
If the particle is close to A scatterers simultaneously before time T, then 

supAA^^"(x) > A. 

Tile r by squares 7f with side length 2ne" and i = 1, . . . , [(l + 4A£"B)T/ne + l]^ 
Then 

{supA/;'"(a;) > A} C {sup sup A/'^^"(x) > A}. 

For each i, cover 7? symmetrically by the larger square 7? of side- length 4n£; 
this implies that 

Hence, by the Poisson field's translation invariance, 

P^({sup > A}) = ((i±i^£!^ + 1)' P^(AA(7f ) > A). 
t<T \ ne J 

Markov's inequality with an exponential finishes the proof. □ 

3.3. Bounding the crossing times. We would like to control, again for a E 
(0, 1/2), the time that the particle takes to cross the ball S^{xf"), compared to a 
trajectory crossing the ball with constant velocity. Define: 

Ux := {supAr„%"i(xr") < A} n {sup|<'"| > 1/2}. 

s<T s<T 

We will call the first set in this intersection Afj_|_]^, and we write B' := sup|F(a;)|. 

Lemma 3.2. Ift'^{s) andtf^^s) are the first entrance and exit times of the trajectory 
with respect to the ball S^{xl), then for trajectories in Ux and for e < {16n\B')~^^°' , 
we have: 

(3.3) il{s)~tl{s)<Ans", 

(3.4) sup sup \vUt) - vUtUm = Ois'^). 
s<rt^(s)<t<t= (s) 

Proof. If a particle enters 5^(a;^'") with constant velocity vl'" = u^(i^(s))/2, posi- 
tion x^{t'^{s)), at time t'^{s), then its exit time i^(s) will satisfy 

iUs)-€{s)<2ne\vl\-' ^ine"^. 
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By the point process description, (3.2), we have for t < t^(s): 

K{t)-vi{tum< E r Ki<{t')-r)dt' 

Putting these two facts together, we have for t < i^(s): 

Kit) - vr\ < M^tMne^'e^'^-'B'Ane- 

So on U\, and for e < (16nAi?')^^/", we have: 



□ 



3.4. Tightness. Wc now turn to showing tightness for D^, the family of measures 
induced by the cut-off processes Uj'" — v^{t At"). 

Lemma 3.3. The family is tight in C{[0,T];M? x M."^): that is, for all 7,77 > 0, 
there exists 6 such that for e small enough, 

pA sup >7l <r?. 

\t — s\<5 

Proof. We use the point process description (3.2) to write: 

(3.5) - ^ A + B + C + D, 
where 

(3.6) ^■■= K{da,du) F^^{x-^%{t'))dt' , 

JSJsAt^ Ju 

(3.7) B:= K{da,du) F^{xi'"At'))dt' , 

J S J S/\T^ Jt'^ ^ 

(3.8) C':=// K{da,du) Fl{xT„{t'))dt' , 

(3.9) D:= Av^'"(i)-Ai;^'"(s). 

Here we have defined t^ g. := inf{t > u : x^^{t) ^ S(a:^'" — cre,£)} A r^, i.e., the 
particle's first exit time after u. Terms B, C, and D concern scatterers that the 
particle encounters at self-crossings, at time s A t^, and at time 0, respectively. By 
the overlap and crossing-time lemmas, 3.f and 3.2, we can bound these terms: 

(3.10) sup\B + C + D\<AKB'nXe°'. 

s,t<T 

If A = e^^, with /3 < 2a, then these terms are negligible in the limit e ^ 0. 

To handle A, on the other hand, requires the martingale decomposition, (3.1). 
However, we run into a problem trying to implement this decomposition, because 
the function f{a,u) := J^"'" F^(x'^'^(t'))dt' , is not adapted to the family of sigma 
algebras, J^t, generated by A'^^-instead, it's anticipative. To replace / by an adapted 
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function, we do a Taylor expansion of around the line i^'" (t) = — u) + ae, 

for t € [u, tfi^ where ^ is the first re-entrance time after u: 

il, := inf{t > u : G - <je,e)} A 

Then the Taylor expansion is: / // 
(3.11) 



F'^{x^%{t'))dt' + / F-^{xt;%{t'))dt' + III + jy, 



where /// is the linear term, and IV is the remainder, with mean value y^{u, a, t'): 



III 



{x'C^t')-xT^{t'))-SIFl{il,^{t'))dt', 



IV 



'di^TAt')iTAt')) ■ yfF^{y^u,a,t'))dt' 



We observe that term / is J^u-adapted and continuous in u. A useful estimate 
from the crossing time lemma and sup„<T.E ^ ^ u\ < Ae is: 

sup <C^A£i+". 

ti<t'<tf,_^ 

The other terms are still anticipativc, so they need to be analyzed further. For 
terms // and IV, it suffices (by sending A ^ oo in the overlap lemma) to show that 
they decay asymptotically on the set 

A^={supAAr(^r)<A}. 

s<T 

In order to bound II and IV, we use the crossing time lemma, 3.2, to show that 
on Afj, both of the following are of order X'^e^'^": 

sup / 

s<t<TJsAT' JS 



sup 

s<t<TJsAT 



Kida,du) 



K{d<j,du) I i^„-«-«(t'))di', 
1 



:(«':(^')i^'.:(i'))-V)^^^„^(2/^(^,a,t'))dt' 



These terms go to zero as £ — > 0, for A = e with /3 < 1/3. Hence, for term IV, 
we can estimate: 



sup 



sAr' JS 



K{da, du) 



i([<-(i') - xf-:(t')] • ^YFl{f{u,cj,t'))dt' 



< 



TAt' 



Nl{da,du) CA^e^" sup|ei • Ve2 • Vi^| 



s Jo 



Choosing A = e ^, with j3 < 3a, implies that term IV is negligible as e ^ 0. 
As for term //, observe that on A!, , for all u. 



(3.12) 



t — u 



-ds 



< CA£i+". 
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From the crossing time lemma, we have — u < Ce, and hence 



sup <C'Ae° 

u<t<K „ 



Then — t^^^\ < C"Ae^+", so we can deduce a version of (3.12), uniformly in t: 
sup <CAei+". 



Then since 



sup 



sup \F^iy)\^B'C"X, 

|y-aS«|<C"Ae2+2a 



< N'^{T)C"^B'X'^s^+°' 



Now to estimate ///, decompose v^'" much like [DGL]: 
(3.13) vr^yl{t)+yt{t)+y^oit)- 
The three processes 2/+ _ q ^^"^ defined as follows: 



(3.14) 



(3.15) 



Voit) 



s Jo 



Nl{da,du) / ' Fl{il(t'))dt' , 



yUt)-y±[^) ■■= I n K[da,du)^vT{u,a). 



Here we abbreviate dt = dt"'dt"dt' and define 



Aw^'"(w, a) 



u J u J u 



p%{da',du') f ^^^(i:',>(i"'))-VF^(£^„,:(t'))rft 

J u' 

•t' nt" 



Nl{da',du')x[tl,>u) 



u' J u J u 



i^.^«'':.'(i"'))-VF„-(C;:(i'))dt 



(3.16) 



t ru 



u J u — Ae-^ J S 



Nl{da',du') / i^^«C'(^"'))-V^^^(i^„,:(i'))rft 



Az;l(M,cr) ■ 



U ^ U 



AN^u) / F^(i^f,,(t'"))-VF^(x^::(t'))dt. 



Noting that Av^°'{u,a) is already adapted and left-continuous, we have further 
split Awl'" into its adapted left-continuous part, Awl'", and a jump part, with 
AN^iu) := N^{S, [0, u]) ~ N^,{S, [0, u)). 
It remains to show that: 

(i) There is a positive constant C such that for every e ^ 1, 

E^(|2/|(i) - y|(s)n < C\t ~ for \t ^ s\ < e; 

(3.17) (ii) For every 7, 77 > and £ <C 1, 

P%{ sup \yUt) - yUs)\ > 1}) < V- 

s,t<T 
\s-t\<e^ 
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Similar to [DGL], the proof uses the martingale splitting developed earlier, as well 
as the Burkholder-Davis-Gundy inequalities. For example, to obtain the first in- 
equality of (3.17) for yl, we use the martingale-compensator splitting: 



E-{\yt{t)-yt{s)\)=E- 



N'Avt + 



1 2 



(3.18) 



< 4 ME^ 



//.;Aa)1+r[(//M;Aa) 



F' ■ VF' 



We also need a bound on the Poisson rate: p%{da, du) < ^p'^^edadu = |p£ ^"dcrdu. 
Then is dominated by N^, which is Poisson with density ^TTpe~'^"dadu. And 
for the first term on the right-hand side of (3.18): 



E^ 



pAvl 



< E 



p)\Avty 



Then use the quadratic variation formula and continuity of p"^ to bound the 
second term by the first: 



E^ 



p^Av^ 



= E' 



p^Avl 



Thus the first two terms of (3.18) arc bounded by a multiple of \t — sp. The 
other terms can be handled similarly, as can the inequality for y^, since 



sup |t^,„ 



u\ < 4e. 



Next we handle tightness for yQ, showing that for some positive constant C, 

E^i\ym-ym\'^'') <C'\t- s\'+''/\ 

Let K be the number of self-crossings in the trajectory, and let be the tube 
around the trajectory up to the stopping time r^. At the z-th self-crossing, let 
Tf and Tf be the times that the particle enters and exits the tube. Set U :~ 
[0,Tf] U [Tf,T|] U ••• U [f^_;^,r|,], and Poi% is a Poisson point measure with 
intensity p%{da,du) := —pe~'^"v^'°' ■ dcrdu V 0. Then define: 



N^ida,du) :-- 



N^{da,du) on [/, 



Poi%{da,du) oSU. 
Then 3.2 implies that there is a number n = n{(j}^ a) independent of e and such that 



s<T 



The remainder of the proof, showing tightness for yQ(t) :— JJ J F^{xf^'^) is a 
straightforward generalization of arguments in [DGL]. □ 
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4. Characterization of the limit 

To identify the limiting process as the velocity diffusion associated to the linear 
Landau equation, we use the Stroock-Varadhan martingale formulation, thereby 
showing that the tight family of processes converges to the velocity diffusion, we 
show that two quantities involving the process and the infinitesimal generator of 
the diffusion are martingales; the desired result will follow by Levy's lemma. 

The infinitesimal generator of the diffusion process is: 

(4.1) L = np\/pJ (fk{k®k)S{k-p)\V{\k\)\'^-\/p. 

It can be expressed more simply in polar coordinates, p = ir,9), as the Laplace- 
Beltrami operator: 

(4.2) L = (:d^/dd^. 

Lemma 4.1. Let s,t S [0,r], f be a smooth test function, and (f>s be a smooth and 
bounded test function that depends only on p{u),u < s, where p G C([0, T]^; M^xM^). 
Also write = (j)s{v^'"). Then we have 



(4.3) lim 







t/\T' \ 



0; 



Proof. It suffices to prove (4.3) for linear and quadratic functions /; we restrict our 
attention here to the linear case. Using the decomposition (3.5) we can see that the 
B, C, and D terms are negligible in the limit, as follows. We cut their expectation 
into two pieces, S := {sups<j- A/'^j^ a)('^) — ''^J' ^^'^ ^^'-^ '^^ ^^^^ term we use 
the bound on the supremum of B + C + D from (3.10), and on the second term we 
use the overlap lemma, (3.1): 

W[iB + C + D)cj,l] < CA£" +C're-iE-[supAA„^(^.,)(,s)x{supACVa)(s) > A}] 

s<T s<T 

< CAe" + C'Te-i / P-{sup AC!(^,,)(s) > X'}dX' 

< CXe^ + C'Te-^ (l±^^^^^ + ' cxp i2n^pe~ 



-2a-l 



X'dX' 



< CXe°' +C'Te~^h{X,e,T)e 



-A 



Here, h involves terms like A^T^e^"^^, and A can be chosen to be a small enough 
(for example, £-"+3/2+/32^-i-/3^ with f3 small and positive) to get the right-hand 
side to vanish with e. Hence only the term A matters: 




K{da,du) / F-{xl;%{t'))dt'\^l 



S J S/\T^ 



Next recall the decomposition of (3.13), Wj'" ^ y+it) + V-it) + J/o(0- We define 
yQ{t) and yj.{t) by replacing by N^, the point process that ignores self-crossings, 
in their formulas; this results in negligible errors. Then we can combine this with 
the Taylor expansion in (3.11), and use the martingale property of to arrive at: 

(4.4) w[{vr - v's'^^Ws] - E'W-it) - ytis) + m) - yu^m]- 
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Then we can use the formulas from the previous section, (3.14) and (3.16), to reduce 
the proof of the lemma to the following two limits: 



(4.5) 



lim 



lim 



'S J S/\T- 

rt/\T' 



N'Jd<j,du)[Adtiu,<j) + Av%{u,a)] j cbl 

ft' nt" 



Nl{da, du) 



S J sAt^ 



u J u J u 



F^{xli%{n)-yFi{xl:'l{t'))\cj^\ 









[( 








J sAt^ J 





lim E^ 



The first of these two limits can be handled by manipulations similar to [DGL], 
except with pl{da\du') = -~pe~'^°'vf^'"' ■ daduW and ij;',%,_„(t) wf;"(t-M') + cr'e. 

For the second of the two limits in (4.5), we examine the part without the infini- 
tesimal operator. Using the martingale decomposition of iV^ and the definitions of 
p, F^, and Vi^^, wc arrive at: 



E^ 



Kida, du) 



S Js/\T^ 



t' rt" 



u J u J u 



Fi{xT.{n)-wFi{ii:%{t'))\cj^: 



S J sAt' 
sAt^ 



p%{da,du) 



t' pt' 



U 'J U 'J u 



^^(i^.::(t"'))-vF„-(x^'.:(0) 



pe-^^du I v^'^da 



it „ rt' rt" 



U 'J U 'J u 



■e°'-^VF 



Now we manipulate the inner integrals to make them look like the desired operator: 
(4.6) 

- / p£"2" / t")Pi^u"t" + a) ■ VF«-"t' + a)dt"dt'dadu 



ps 



'2a 



2^' 



-2a 



TF{r + ijf;"r) • \7 F{r)dTd^rdu 

V? • F{r + Cr)F(r)] drd^rdu 

e''=-^^(fc(8)fc)|F(|fc|)p drd^kdu 



sAt' 
tAT^ 



sAt^ 



p£~^"V^ • / TTS{k-Oik(g>k)\Vi\k\)\^ d^kdu 



e-^^Lvt'^du. 



□ 
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5. Convergence in expectation 

We will show that, in particular, for an initial distribution /o, 

(5.1) /^(t, X, v) = E[U^l^Jx, v))] ^ h{t, X, v), 

where h is the solution of the linear Landau equation, (1.3). Although this is a 
weaker mode of convergence than that just proved, the following argument illus- 
trates the key intuition behind the result for the entire range of a: Although there 
are significant numbers of clusters of obstacles, their total influence is actually neg- 
ligible. The initial set-up is similar to [DR.] but differs after the first two steps. 

First, observe that f^{t,x,v) can be written as (using |-B(a;)| to denote the mea- 
sure of the ball of radius t around the initial position x): 

/^(t,x,z;)=e-''^l^(-)I^Mli! / ... / M<^l^Jx,v))dn...dr^. 



uj.e \ 

B{x) JB(x) 



Define a cut-off xi j killing configurations that have an obstacle at the initial position: 

Xi{u) :=x({^ = {'^aili :Vi = l,...,iV,|x-r,| >e}). 

Making this cut-off introduces an asymptotically vanishing error. That is, there 
exists a function 0i(£) as e such that can be written as: 

Next we define another cut-off, X2^ killing configurations with obstacles that arc not 
encountered by the light particle and thus have no effect on the tra,jectory: 

X2{oj) := x({c^ = {n}-U : V* = 1, . . . , Af, r, e T{t)}). 

Here T{t) is the tube of radius e around the light particle's trajectory: 

T{t) :=={2/:3se [0, i] such that \y ~ x{s)\ < e}. 

Then we have again an asymptotically vanishing error: 

Next we make the key observation that single obstacles dominate the trajectory's 
path. First, we note that there is a sequence of thresholds, q;„ converging to 1/2 
from below, such that for a < «„, there is a negligible number of clusters of n 
obstacles. For instance, a2 = 1/4, that is, only for a > 1/4 do we need to worry 
about the numbers (if not influence) of doublets (clusters of 2 obstacles). 

Let n be the number of internal doublets up to time t; n is a random variable 
with an expected value of order e^~^". Let 6j be the deflection angle for the collision 
with the j-th doublet. 

Then the key observation is that E'^[X]j=i ^j] ^ as e ^ 0, and can be seen 
as follows. We introduce an expansion of the deflection angle at the j-th doublet 
(which comes from (6.7) in the next section): 

=£"A{+e2"A^2+0(£3"). 

We note that E0j = 0, even when the test particle hits a cluster of rij scatterers 
overlapping to form one large obstacle. This can be seen easily for the cluster 
which is actually just one scatterer, because 9j is an odd function of the impact 



/^(t,a:,z;) =e-^=l^(*)l ^ ^^TT / XiX2(t^)/o(<I>Ue(^, + 02(e) 
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parameter. For nj > 1, we assume without loss of generality that after leaving the 
previous cluster, the light particle's position is xj = 0, and its velocity, Vj = ei 
(the first coordinate vector). Let ri be the centers of the scatterers in the order 
that they are encountered, with i between 1 and nj. For each such configuration 
Rj := {ri}^jLi, there is an opposite configuration Rj, obtained by reflecting across 
the first coordinate axis, with corresponding scattering angle 6j = ^Oj- So 9j is an 
odd function in this sense, and vanishes when integrated against an even measure. 
Then we compute: 



E 



E 



E[n]E[e2"A^2] 
0(ei-^")0(e"") 



And this last quantity goes to zero as e because a £ (0, 1/2). 




Figure 4. The changes of variables, from r.j to tj and bj, and from 
hj to 9j. 



Resuming along the lines of [DR], we introduce a change of variables £, replacing 
obstacle locations by hitting times and impact parameters; 

(5.2) C:{T,Y^^,^{t,,b,}^^,. 

The difference is that here, we define the domain of C to be the set of all trajectories 
that do not start on a scattercr; that do not have extraneous scatterers; and that 
have stopping time > t (i.e., they have no almost-tangential self-intersections, 
not too many self- intersections, and velocity that does not change too much). In 
particular, we handle the change of variables at a doublet as illustrated in figure 5. 
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Figure 5. The change of variables for a doublet. 



Then we can write using the change of variables and (2.7), with an error 03(e) 
vanishing as e — > 0, as: 



M>0 ■'^ 



Here t stands in for (ti,...,tfc). and A {t ; ti G (0, i),t2 G {ti,t), . . . ,tk G 
itk-i,t)}. Similarly, b := {hi, . . . ,5^), and (-£,£) x • • • x (-£,£). 

Then the procedure of [DR] can be followed and a change of variables made (with 
a slight modification for doublets, described below), from impact parameters {bi} 
to deflection angles {9i} (see figure 4), which has Jacobian determinant: 



M M 

n.-"r.(..):^nt- 

i=l 1=1 * 



Here, ^/.{Oi) is the rescaled scattering cross section. We use the fact that the 
deflection angle 9 through a doublet can be approximated: 9 = 9i + 62 + (/»(£)■ Here 
9i and 92 are the deflection angles corresponding to bi and 62 in figure 5, and 
is a small error vanishing as e ^ 0. 

Using also i?e(w) to denote the rotation of the vector v by angle 9, we make the 
following polygonal approximation to the trajectory: 



M 

x{t) = X + J2Re,+-+0Av)it'^+i - U) + 0{Me). 

i=0 
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Additionally, we approximate by 2et. Then we can rewrite 

X /o (^x + ^ Rg^+...+0^(v){U+i - ti), d9dt + 04(e) 

M 



j=0 ^ 



M>0 "A^n^j^j 



/ ^ \ - _ 

X /o ( a; + ^ Roi+.-.+e, {v){tt+i - U), Rei+---+eM (^^) ] dOdt + 04(e) 

We then recognize this expansion as the series form of a solution to the family of 
Boltzmann equations: 

(dt+v ■V^)h,{t,x,v) ^ p I T,{e)[h,{t,x,Re{v)) - hS.x,v)]de 

J —TT 

he{0,x,v) = fo{x,v). 

Finally, the converge in the appropriate sense to h, the solution of the Lan- 
dau equation (1.3), because the scattering cross sections concentrate on grazing 
collisions. 

We note that E6j — 0, even when the test particle hits a cluster of nj scatterers 
overlapping to form one large obstacle. This can be seen easily for the case of 
hitting a lone scatterer, because 9j is an odd function of the impact parameter. 
For Tij > 1, we assume without loss of generality that after leaving the previous 
cluster of scatterers, the light particle's position is Xj ~ 0, and its velocity, Vj = ei 
(the first coordinate vector). Let r, be the centers of the scatterers in the order 
that they are encountered, with i between 1 and nj. For each such configuration 
Rj there is an opposite configuration Rj, obtained by reflecting across 

the first coordinate axis, with corresponding scattering angle 9j ~ —Oj- So 6j is an 
odd function in this sense, and when integrated against an even measure, vanishes. 

6. The diffusion constant 
We take the following as the definition of the diffusion constant: 

9^r^{0)de ^ \ime-^"- e{hfdh. 

-TT 2 

Proposition 2. The diffusion constant, Q, defined above, is independent of a & 
(0, 1/2) and can he expressed by the following formula: 



2 7-1 \Jb \u J u vT^ 



u 



Proof. We do several expansions to compute the diffusion constant. First, we com- 
pare the deflection through one obstacle, centered at the origin, to the path that 
would be taken if the obstacle weren't there. Let the entry time, place, and velocity 
be 0, x~ = a;(0), and v~ ] the deflected exit time, position, and velocity be t, x'^ , 
v'^] and the non-deflected (straight-line) exit time, position, and velocity be f, i"*", 
and . (See figure 6.) 
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Figure 6. The straight-line trajectory exiting at time f and the 
deflected trajectory, at r, with the obstacle's center at the origin. 



We desire an estimate for |r — fj. so we examine the identity: 
£2 = |^+|2 ^ \x-\'^ +t'^ + 2tx- -y- 
= £^ + + 2tx~ ■ v~ . 

This identity gives f = —2x^ ■ = 0(e), as we may assume that \v'\ — 1. 
Similarly, r = —2x^ ■ {v^ + 0(e")), so we have the estimate: 

(6.3) \r-r\=0{e'+^'). 

Next, we estimate u+ — v~ , using the definition F^{y) :— ~W^{y). We could 
make the following estimate (but wc will actually do better): 

F^{x{s))ds 

f F^{x{s))ds+ [ F^(x{s))ds 
Jf 

F^{x{s))ds + 0{e^+'') ■ 0{e"-^) 

F^{x{s))ds + 0{e^''). 

Accordingly, we estimate that for s between r and f, and some a between r and s, 
and using the compact support of V: 

F^i^is)) = F^{x{t)) + VF^(a)(x(s) - x{t)) 

(6.4) =0 + O(£"-2) •0(ei+") 

= 0(£-i+2"). 

We now combine (6.4) with (6.3), the estimate for |t — t|. Instead of the previous 
remainder term, 0(6^"), we get: 

F^{xis))ds ^ 0(£i+") • 0(£"2+") . 0{e^+") = 0(£3"). 

J f 

Thus: 

(6.5) v+-v-^f F^(x(s))ds + 0(£3"). 



18 



KAY KIRKPATRICK 



Expand F^(x(s)), using x{s) = x + tv , and rename the first two terms: 
= KiHs)) + VF^ixis))ix{s) - xis)) 

+ ^D^F^ia)ixis)-xis))-{xis)-x{s)) 
=■■ Y1+Y2 + Ci(e-3+")Ci(ei+")2 
To analyze Y2 further, observe that: 

yF^{x{s)){x{s)-x{s))=^K{m} I Ht)-v-]dt 

Jo 

= \7F^{x{s)) r I F^{x{u))dudt 
Jo Jo 

+ VF-(x(s)) /' [\fI{x{u)) - F-^{x{u))]dudt 
Jo Jo 

= : Y21 + Y22- 
Regarding the first term, ^21, we integrate out t: 

X21 := / Y2ids = f VF^{x{s)) f f F'^{x{u))dudtds 
Jo Jo Jo Jo 

(s - u)WF^{x{s))F^{x{u))duds 

Jo 

= 0{e^)0{e°'-^)0{e°'-^) = ©(e^"). 
As for the second term in this decomposition, 1^22: 



Y22 = VFl{x{s)) r ( F^ixiu)) - F^ixiu))d 
Jo Jo 

= 0(e~'+")0(e2)0(£ 
Then integrating, we get: 



u 



X22 ■■= I Y22ds = 0(e3"). 
Jo 



Thus we have, combining the estimates for Yi, 121, and I22: 

Xi+X2i + 0{e^"). 




Figure 7. Change of variables from s to which is the distance 
remaining to the center of the crossing (scaled by e). 
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Nowwe claim that := Jq Yi = /J" J^^(i(s))ds is perpendicular to Without 
loss of generality we may assume that v~ = (1,0). Making the following change of 
variables, (s, b) i— > (^, 6), as in figure 7, we see that the claim is true: 



Yi = - I V 







\x + sv \ \ X + sv 



\x + sv \ 

1 



ds 



^ /_J'(v/5^)--^(^,^)rfe 
Then we introduce a variable u such that: 

M = and thus ^du=^mL. 



Using this change of variables, and uq := \b\/ ^if- ^IP- = \b\ 

1 



















J -1) 



(6.6) 



2e" / V'{^¥+e) 



1 



bd^ 



2e° 



,f\b\ 



V'{'-^\^b{- 



uj \b\ 



\b\' 



\b\VT 



du 



„ / , f\b\\b du 



u J u - V? 



So if wc examine the deflection angle, 0, as pictured in figure 6, we obtain: 

1 - + 0(e^") = cos6i = t;+ • = {v+ - v^) • + 1. 
Then orthogonality of Xi with implies: 

= -2 /" • VF^(i(s)) ( [s ~ u)F^{x{u))duds 
Jo Jo 

= -2^ 1f;(£(s))^ {s - u)F^{x{u))duds 

= -2 F^{x{s)) r {s ~ u)F^{x{iL))du 
Jo 

d 



(6.7) 



+ / ^a(i(s))^ / (s - u)F^ixiu))duds 



= + 2/ F^(x(s)) / i^^(i(u))(iMds 

= 2 / i^^(x(s))i^^(x(M))x{„<4(u,s)dwds 
Jo 

2 

F^(x(s))ds 



= 1^1^ 
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We note that in particular, (6.7) gives us an expansion of the deflection angle with 
leading term of order e". In terms of the difference in velocities: 

v+ -V- =Xi +X21 +0(e^") 
1, 



(6.8) 



-Xi-||Xi|2 + 0(e3") 

= Xi- ^6(2 +0(£3") 

= 0(e") + 0(e2") + 0(£='"). 



The proof is finished by putting this expression for 9^ into the definition of the 
diffusion constant, (6.1), and using (6.6). □ 

Lemma 6.1. The diffusion constant for a G (0, 1/2) is the same as for the case 
a~\/2 in [DGL], whose formula via the martingale characterization was: 

(6.9) C:^np\vm-' l\kf\V{\k\)fd\k\. 

Proof. This can been seen by comparing the previous computation with (4.6). We 
take the quantity J (Yi + 121) and integrate it with respect to x~ , change variables 
{t = u — s and r = + w^s), and apply Plancherel's theorem: 



{Yi+Y2i)dsdx =11 I {s-u)F^{x +uv^)VF^{x + sv^)dudsdx~ 
J Jq Jo 

= J J TF^{r + TV-)\/F^{r)dTd'^r 

= \j j [{'^p-F{r+pT))F{r)dTd\]^^^_ 



Vp I I e'^-P{k®k)dTd^k tp 



7. Appendix 



□ 



Retaining the same essential argument, the result of [DR] (convergence in ex- 
pectation of the evolution of the initial distribution /g) can be stretched to include 
a G [1/8,1/4). The estimate that must be improved in order to achieve this is 
(40) in [DR, Lemma 1]. This section illustrates how to iterate their geometric 
method to get a tighter bound on J"^ , the error term estimating the probability of 
non-consecutive overlappings and recollisions, a bound that decays for a < 1/4. 

First fix n disjoint closed subintervals of (0, tt), each of the form 

Im [4>rn, 4>m + 7r/2n]. 

We assume that < (f)2 < ■ ■ ■ < (t>n- For this iterative method to work, n must 
satisfy 7r/2ri > Ce" > |6'fc|. Thus for each m = 1, ... ,71, there exist {/i,ri}"i=i such 
that 

0k G Im- 

k=l 

Our trajectory will have n gaps where the times th„^ can vary. All other times 
tk and all angles 6k arc fixed. 

The case n = 2, /i = [tt/S, 37r/8], and I2 = [57r/8, 77r/8] is illustrated in figure 8. 
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Figure 8. A recoUision {x{tj) lying in the tube of radius e around 
the trajectory) constrains how much hi and /12 can vary in total; 
the variation for each is bounded by \/2e. 



We claim that the recollision condition. Vj e Use(t t 1) -^(^(■5)1 2e); together 
with each Im being bounded away from the endpoints of [0,7r], results in each t^^ 
taking values in a set whose measure is 0{e). The claim can be seen by working 
backwards from /i„: the height of the trajectory at time t^^ can vary only in an 
interval of size e, due to the recollision condition restricting the trajectory to a tube 
of width £. Then the total variation of all the t^^ , m = 1, . . . ,n can be no larger 
than Ce. Here the constant C is chosen so that 



sin( 



< C, for e {0i,0„ +7r/2ri}. 



Hence the variation of any single is bounded above by Ce, there being no 
negative values of th,„ allowed and all angles pointing upwards at the variable times 
(i.e., being strictly between zero and tt). Then J"^ can be estimated like before: 
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1 Q 



Q>1 "'tQ-l 1=0 j=!:+2/ii=i+l 

E ■•■ E ni((E'^'^-e^"|)'((^^^ U ^(c(-)^2e)y 

/t2=i+2 h„=i+nm=l fc=l sG(ti,ti+i) J' 



Q>1 

< (^^y^£3n+2 — (2n+2)(5 ^^jn-jgn— 4(n+l)a 

The exponent is positive if a < , so taking n to infinity (as e goes to zero) 

gives the resuh for a < 1/4. 
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